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Abstract . Electronic transport in a graphene-based ferromagnetic /normal / ferromagnetic 
junction is investigated by means of Landauer-Biittiker formulism and the nonequilib- 
rium Green's function technique. For the zigzag edge case, the results show that the 
conductance is always larger than e 2 /h for the parallel configuration of lead magneti- 
zations, but for the antiparallel configuration the conductance becomes zero because of 
the band-selective rule. So a magnetoresistance (MR) plateau emerges with the value 
100% when the Fermi energy is located around the Dirac point. Besides, choosing nar- 
rower graphene ribbons can obtain the wider 100% MR plateaus and the length change 
of the central graphene region does not affect the 100% MR plateaus. Although the 
disorder will reduce the MR plateau, the plateau value can be still kept about 50% 
even in a large disorder strength case. In addition, when the magnetizations of the left 
and right leads have a relative angle, the conductance changes as a cosine function of 
the angle. What is more, for the armchair edge case, the MR is usually small. So, it is 
more favorable to fabricate the graphene-based spin valve device by using the zigzag 
edge graphene ribbon. 
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1. Introduction 

Graphene, a two-dimensional single-layer crystal of carbon atoms arraying in a 
honeycomb lattice, is a novel and exciting material pp. The low energy excitation 
of a graphene has a linear dispersion relation, and the dynamics of the charge carriers 
obeys the massless Dirac-like equation [2]. Because of the unique band structures, 
the graphene has many peculiar properties, such as that the Hall plateaus occur at 
half-integer multiples of ge 2 /h with the spin and valley degeneracy g = 4, and the 
conductivity at the zero magnetic field has a non-zero minimal value [3j . For the neutral 
graphene, its Fermi level is located at the Dirac points, the corners of the hexagonal first 
Brillouin zone. In experiment, a gate voltage can be used to tune the Fermi level, which 
can be above or below the Dirac points [I]. In addition, graphene also exhibits many 
excellent transport characteristics: the high mobility [5j, long spin relaxation length [6], 
and stable behaviors under ambient conditions. At room temperature, its mobilities can 
be above 10 4 cm 2 /V • s, implying that the mean free path can be as long as a few hundred 
nanometers. Because of weak spin-orbit coupling [7] and low hyperfine interaction [8], its 
spin relaxation length can reach the order of the micrometer at room temperature. Thus, 
graphene may be an excellent candidate for microelectronic applications, in particular 
for spintronic applications [9]. 

In general, the charge carriers in graphene are not spin-polarized. For spintronic 
applications, people try to inject the spin-polarized current or to induce the spin- 
polarized carriers in the graphene. Recently, many works are focusing on this issue. For 
example, Haugen et al |10| suggested that the spin-polarized carriers can be realized 
by growing the graphene on a ferromagnetic (FM) insulator (e.g. EuO). Owing to the 
magnetic proximity effect, an exchange split between the spin-up and spin-down carriers 
in the graphene is induced, and then the carriers are spin-polarized. Based on the first- 
principles calculations, Son et al [11 J predicted that the zigzag graphene nanoribbon 
becomes a half metal when an in-plane transverse electric field is applied. Also, Lin 
et al [12] demonstrated that electron-electron correlation in armchair graphene nano- 
ribbon can generate the flat-band ferromagnetism. On the experiment side, a large 
spin injection into the graphene has been realized by connecting the graphene to an 
FM electrode [El [131 [H]. Furthermore, several groups [HI [15] have performed nonlocal 
magnetoresistance(MR) measurements, in which a net spin current is brought between 
the injector and detector. 

The well-known application in spintronics is the spin- valve effect [16], in which 
the resistance of devices can be changed by manipulating the relative orientation of 
the magnetizations. Motivated by the spintronic application with the novel material, 
the spin-polarized transport through graphene are currently attracting a great deal of 
attentions pU EU CEH1 EDI ED E21 E3] . Using the tight-binding model, Brey and Fertig 
[20] studied the MR of the FM/graphene/FM junction in the limit of infinite width. 
They found that the MR is rather small since the conductivity is weakly dependent on 
the relative magnetization orientations of FM leads. Ding et al [21] investigated the 



Electronic transport through a graphene-based ferromagnetic/normal/ 'ferromagnetic junction?) 

similar device with two FM leads by a continuous model. The results showed that the 
MR versus the bias exhibits a cusp around the zero bias in absence of external magnetic 
field and oscillating behavior at the high magnetic field. Based on the first-principles 
calculation, Kim and Kim [22J predicted that a graphene-based spin-valve device could 
have a high MR. 

In this paper, we study the conductance and MR of a graphene-based spin valve. 
The device consists of a graphene nanoribbon coupling to two FM leads, as shown in 
figure [U^a). Here the width of the device is finite in the order of 10 nanometers, i.e. 
the size effect is considered. In recent experiments [23], a few 10-nanometer or sub- 
10-nanometer graphene nanoribbons have already been fabricated successfully. For a 
finite-width graphene nanoribbon, the wave vectors along the confined direction are 
discrete, and the transverse subbands emerge. The characteristics of the subbands are 
strongly dependent on the chirality of the graphene-nanoribbon edge |25j, e.g. the zigzag 
edge or armchair edge. So the conductance and MR should be strongly dependent on 
the boundary condition of graphene nanoribbon. In addition, we consider that two FM 
leads also have the hexagonal lattice structure as the graphene. In other words, the 
leads are graphene-based FM or called the FM graphene. On the experiment, the FM 
electrode (e.g. the cobalt electrode) usually overlaps on the graphene through a thin 
oxide layer [6j[Tl]. Because of the magnetic proximity effect and the Zeeman effect, an 
exchange split in the graphene underneath the oxide layer is induced, then its carriers are 
spin-polarized. So the graphene covered with the FM electrode has the magnetization. 
Then the spin-polarized charge carriers, driven by a bias, travel from one FM graphene 
through the central normal graphene to another FM graphene. 

In the tight-binding model the Landauer-Biittiker formula and non-equilibrium 
Green's function method are applied to calculate the conductance and MR. For the 
zigzag graphene ribbon case, we found that when the graphene nanoribbon is narrow 
enough and its separation A between the first subband and the zeroth subband is larger 
than the exchange split energy M, the conductance for the antiparallel magnetization 
configuration can almost be zero in a quite large energy window around the Dirac point. 
But for the parallel configuration the conductance is always larger than e 2 /h, so the MR 
exhibits a plateau with the value 100%. As the width of the nanoribbon increases, the 
sub-band separation A gradually decreases, and the MR can keep at the value 100% at 
the beginning of A > M, then decreases when A < M. In the presence of the disorder, 
the MR slightly drops, but even in quite large disorder it can keep the value 50%, which 
is still much larger than 10%, the lower limit of the giant magnetoresistance effect. What 
is more, for the armchair graphene nanoribbon case, the MR is always small regardless 
of the width of nanoribbon. By comparing two cases of graphene ribbons with different 
edges, it is more reasonable to apply zigzag graphene ribbons to the spin-valve devices. 

The rest of the paper is organized as follows: In Section |2} we describe the model 
and the details of the calculations. In Sections [3] and H] we will give the numerical results 
of the conductance and MR for the zigzag and armchair edge cases, respectively. Finally, 
a brief conclusion is presented in Section |5j 




Figure 1. (Color online) (a). The schematic of a graphene-based FM/normal/FM 
junction, (b) and (c): the conductance G vs. the Fermi energy Ep for the different 
magnetization M at the parallel configuration (8 — 0). The size of the central graphene 
is N = 50 and L = 10. 



2. Model and formulations 

We consider a graphene-based spin valve (as shown in figure Ufa)) which consists of a 
central normal graphene strip and two FM graphene ribbons. The total Hamiltonian 
H of the device can be divided into four terms, H = He + Hl + Hr + Ht, where 
Hq describes the central normal graphene region, Hl and Hr are the Hamiltonian 
of the left and right FM-graphene leads, respectively, and Ht = Htl + Htr is the 
coupling Hamiltonian of central region to the left and right leads. In the tight-binding 
approaximationpS] . the Hamiltonians He, Hl, Hr, and Ht can be written as 

He =^2eca\a I a i - (tajaiaj + h.c), 

i£C <ij>(i,jeC) 

H a=L ,R = a l( e a°"/ + CT • M a )aj 

— (ta\aiaj + h.c), 

<ij>(i,jEct) 
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H T = — {ta\aiaj + h.c), 

<ij>(ieC,jeL,R) 

where s =t, 4 represents the spin of electrons; a\ s (a, s ) creates (annihilates) an electron 
with spin s on site i, and Oj = (" i1 "J; <U> stands for a nearest-neighbor pair, 
cr = (o^, o~ y , a z ) are the Pauli matrices and 07 is a 2 x 2 unit matrix, e^, e^, and 
e# are the on-site energies (i.e. the energy of the Dirac point) in the center region, left 
and right FM leads, respectively, which can be tuned by the gate voltage. The size of the 
center graphene region is described by the width N and length L. In figure [T^a), a zigzag 
edge graphene nanoribbon with N = 4 and L = 9 is shown. The terms including the 
factor t in Hamiltonian describe the nearest neighbor hopping with the hopping energy 
t. and are the magnetizations of the left and right FM leads. Here we allow 
that the magnetizations and M# can be along arbitrary direction. Without loss of 
generality we assume that the magnetization of left lead is along the z-axis, then 
M L = Ml(0, 0, 1), and the magnetization M R of the right lead is along the direction 
(9, ip) with = Afft(sin 9 cos (p, sin 9 sin <p, cos 9). 

Before performing the calculation, we take a unitary transformation with 5j = t/Oj 
for all site % in the right FM lead, where the unitary matrix U is 

( cos(fl/2) e-^sin(0/2) \ 
= ^ sin(0/2) - cos(#/2) ) ' 

Under this unitary transformation, the Hamiltonians He, Hl, and Htl remain 
unchanged, and H R and H TR vary into 

— (ta\ajaj + h.c), 

<ij>(i,jeR) 

Htr = — X! (ta\Udj + h.c), 

<ij>(i£C,j€R) 

where M' R = M R (0, 0, 1). After the unitary transformation, the z-axis of the spin in the 
right FM lead is along the direction of M^, and the Hamiltonian H R is diagonal in the 
spin space. 

The current flowing through the device can be calculated from the Landauer- 
Buttiker formula [27] 

I=(e/h) J deT LR (e)[f L (e) - f R (e)}, 

where /l/_r(c) = l/{exp[(e — Hl/r) /ksT] + 1} is the Dirac-Fermi distribution function 
of the left and right FM leads and T LR (e) = Tr[r L G I T R G a ] is the transmission 
coefficient, with the line-width function r a (e) = z[EJj,(e) — S*(e)] and the Green's 
functions G r (e) = [G a (e)]t = l/[e - H c - - S^]. S^ /R is the retarded self-energy 
function coupling to the leads, which has to be calculated numerically by solving the 
surface Green's function of the leads [281 129] • After solving the current /, the linear 
conductance G can be obtained straightforwardly, G = liniy^o dl/dV, with the bias 
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V = hl — fiR. At zero temperature G = (e 2 /h)T(E F ). In the numerical calculation, we 
take the nearest-neighbor hopping energy t — 1 as the energy unit and adopt the Dirac 
point energies = €r = ec = 0. The magnetizations M a in the left and right FM leads 
are assumed to be equal (M L = M R = M), which is usually reasonable when two FM 
leads are made of the same material. The angle <p = because the conductance and 
MR are independent of (p. 



Left Center Right 




Figure 2. (Color online) The energy band structures of the left lead, center region, 
and right lead in the case of the zigzag edge at the parallel configuration (9 — 0). (a) 
< M < A/2 and (b) A/2 < M < A. 
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Figure 3. (Color online) The conductance G vs. the Fermi energy Ep for the 
different M at the antiparallel configuration (6 — tt). The size of system is the same 
as figure [TJ^b) . 
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3. The case of zigzag edge 

In this section, we focus on the zigzag edge graphene-based spin valves. We first 
investigate the conductance G on the parallel and antiparallel configurations, and then 
the MR and the conductance on the arbitrary angle between magnetizations Mi and 

Mr. 

3.1. The conductance for the parallel configuration 

Figures HJb) and (c) show the conductance G versus the Fermi energy Ep for different 
magnetization M of the parallel configuration (i.e. 6 = 0). When the magnetization 
M = 0, the two leads are normal and the whole device is a flawless graphene ribbon. 
In this case, the conductance G exhibits the plateau structure with the plateau values 
at 2,6,10,... (in the unit e 2 /h), i.e. at the half-integer position, g(n + l/2)e 2 /h with 
the degeneracy g = 4, due to the transverse sub-band structures. The step height of 
the value 4e 2 /h results from spin and valley degeneracy and the plateaus at half-integer 
value originates from the fact that the zeroth subband has the only spin degeneracy. 

When the two leads are FM with a non-zero M, the conductance G does not have 
the perfect plateau structure, but still obeys the electron-hole symmetry, G(—Ep) = 
G(Ep) (see figures [T^b) and (c)). In the following, we discuss the conductance G in 
details for < M < A/2 and A/2 < M < A, where A is the energy splitting between 
the first subband and the zeroth subband (see figrue [2]) and A f» 0.045 for N = 50. 
For < M < A/2, the curve of G versus Ep can be partitioned into several regions: 
(0, M), (M, A - M), (A - M, A + M), etc. In the first interval of Ep G (0, M), the 
conductance G oscillates with the value between e 2 /h and 2e 2 /h. In the interval of 
Ep G (M, A — M), the conductance G is exactly equal to 2e 2 /h and forms a plateau. 
In the interval of E F G (A — M, A + M), the conductance approaches the value 4e 2 /h. 
With the rising of Ep further, the conductance becomes larger. 

In order to understand above results, we present the energy bands of the left lead, 
center region, and right lead in figure |2j and the band index is specified in each region. 
Owing to the symmetry of bands, only the part of the moment k > is shown. The 0-th 
and —1-th subbands in the FM leads are nondegenerate, but other subbands are two-fold 
degenerated. Because of the parity conservation of the transverse wave function, the 
electrons belonging to the even (odd) parity subbands in the left FM lead are transported 
only into the even (odd) parity subbands of the right lead, which was demonstrated in a 
very recent letter [30] and it is named the band-selective phenomenon. For convenience, 
we use the symbol k) to denote the channel through which electrons are transported 
from the z-th band of the left lead, through the j-th band in the center region, to the 
fc-th band of the right lead, and T it j^ to denote its transmission coefficient. 

With the aid of energy bands in figure Ufa), we now explain the conductance G in 
figure Q][b) with < M < A/2, (i) When Ep G (0,M), the channel (0,0,0) is open for 
spin-up electrons, and T 0i0i0 is exactly 1 due to lack of the scattering. For spin-down 
electrons, the channel (—1,0,-1) is open, in which the parity of the transverse wave 
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Figure 4. (Color online) The energy band structures of the left lead, center region, 
and right lead for the zigzag edge case at the antiparallel configuration (9 = it), (a) 
< M < A/2 and (b) A/2 < M < A. 

function in the center region mismatches with that of the left and right lead, so the 
scattering exists and < TLi j0l _i < 1. As a result, total conductance oscillates with its 
value between e 2 /h and 2e 2 /h. (ii) When E F e (M, A - M), the channels (0, 0, 0) are 
available for both the spin-up and spin-down electrons, and the conductance G exactly 
equals to 2e 2 /h. (iii) Furthermore, Ep increases to the range of (A — M, A + M). 
Besides the channels (0,0,0) of both spin electrons, the channels (l,j, 1) (J = 0,1) of 
the spin- up take part in the transport, which make the conductance G approximate 
4e 2 /h. 

Now we discuss the conductance in the case of A/2 < M < A. Figure (T^c) 
shows the conductance, and the corresponding band structure is shown in figure Hfb). 
Apparently, the conductance G varies more complex than that inO < M < A/2 
(see figures LTJb) and (c)). Here we are mainly concerned about two energy regions: 
Ep G (0, A — M) and (A — M, M) . When E F 6 (0, A - M) the transport involves the 
channel (0, 0, 0) of spin-up electrons and the channel (—1, 0, —1) of spin-down electrons. 
So the conductance G oscillates between e 2 /h and 2e 2 /h, which is identical with case 
(i) when < M < A/2. For E F e (A - M,M), another channel (1,0,1) of spin-up 
electrons joins into the transport, so the conductance is obviously enhanced and the 
value lies between 2e 2 /h and 4e 2 jh. Further, for Ep > M, because more channels will 
be opened, the conductance usually becomes larger than 4e 2 /h. 

3.2. The conductance for the antiparallel configuration 

Next, we study the antiparallel configuration with 9 = it. The results for the 
conductance G are shown in figure [3j One of the main characteristics is a zero 
conductance when Ep is near Dirac point (i.e. Ep = 0). This is very different from 
the parallel configuration, in which G is always larger than e 2 /h. Let us analyze the 
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conductance in details in the following. Because of G(—E F ) = G(E F ), we only discuss 
the case Ep > 0. (i) For < M < A/2, the conductance G is exactly zero in the range 
E F G (0,M) and becomes 2e 2 /h when Ep G (M, M + A), (ii) For A/2 < M < A, 
G = when Ep G (0, A — M), while G is between and 2e 2 /h when Ep G (A — M, M) 
and is exactly 2e 2 /h when i?F G (M, M + A), (iii) If M increases sequentially (M > A), 
around the Dirac point a small salient appears instead of the zero conductance range. 

These characteristics of conductance can be well understood from their band 
structures in figure HI From the band structures of the antiparallel configuration the 
spin-up and spin-down electrons contribute equally to the conductance, and we only 
discuss the spin-up electrons in the following, (i) We discuss the case ofO<M<A/2, 
and the corresponding energy band structure is illustrated in figure HJ^a). When 
Ep G (0, M), only the channel (0,0,-1) is involved, in which the band is even 
parity in the left lead and is odd in the right lead. Because of the band-selective 
rule, its transmission coefficient 7oo,-i = 0, and G — 0. On the other hand, when 
Ep G (M, A + M), the channel (0,0,0) is available so that G = 2e 2 /h. If further 
increasing Ep, more channels are opened and the conductance is larger, (ii) We analyze 
the case of A/2 < M < A, and the corresponding energy band structure is illustrated 
in figure @Jb). For Ep G (0, A - M), only the channel (0, 0, -1) is involved. So G = 
because of the band-selective rule. When Ep G (A — M,M), the channel (1,0, —1) is 
opened. The parity of the transverse wave function in the center region is different from 
those in the left and right leads. The incident electrons may be scattered and ?i o _i lies 
between and 1. As a result, < G < 2e 2 /h, including the spin-down electrons. By 
increasing Ep to the range (M, A + M), the channel (0,0,0) is opened, and the total 
conductance G = 2e 2 /h. Notice that when Ep is close to but less than A + M, although 
the extra channel (2, 0, 0) is involved, the conductance G still is 2e 2 /h because only the 
subband is available in the right FM lead. 

3.3. The magnetoresistance 

After obtaining the conductances Gp and Ga of parallel and antiparallel configurations, 
the MR, defined as MR = (Gp — Ga)/Gp, can be obtained straightforwardly. Figured] 
shows the MR versus the Fermi energy Ep for different magnetization M. Here the MR 
can be very large as far as 100% when Ep is around the Dirac point. In particular, a 
plateau with the value 100% is clearly exhibited in the curve of MR-Ep. This 100% MR 
plateau origins from the zero conductance in antiparallel configuration (see figure |3]). It is 
noticed that the conductance for the parallel configuration is still quite large (> e 2 /h). 
This means that the present device not only has the large MR but also has a large 
variance of the conductance between the parallel and antiparallel configurations. Thus 
it is a good candidate for the spin valve devices. The width of the 100% MR plateau 
is determined by magnetization M and the energy splitting A between the subbands. 
When < M < A/2, its width is 2M so it increases with the enhancement of M. At 
M = A/2, the plateau width reaches its widest value A. By further increasing M, 
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Figure 5. (Color online) The MR as a function of Ep for different magnetization M. 
The size of the central region is N = 50 and L = 10. 



the plateau width becomes narrow and it is equal to 2(A — M) when A/2 < M < A. 
At M = A, the 100% MR plateau disappears. When M > A a valley instead of the 
plateau appears in the vicinity of Ep = 0. Therefore for a larger A, it is more favorable 
for the formation of the 100% MR plateau. In experiment, the graphene nanoribbon 
with its width of several-ten or sub-ten nanometers has been fabricated [21]. Taking 
a 20nm- width ribbon for an example, its A is about 0.12eV, which is usually much 
larger than M. In addition, when Ep deviates far from (i.e. the region of the 100% 
plateau), the MR may be quite small (see figure [5]). This implies that the capability of 
the graphene-base spin valve is optimal when it works near the Dirac point. 

We study the influence of the size of the ribbon between two FM leads on the 
MR. Figure |6]^a) (Figure E^b)) shows the MR as a function of Ep for different width 
iV (length L) and a fixed L (N). With the increasing of width N, the separation A 
of the transverse subbands decreases monotonously (with A w Snt/AN) [31] . In the 
beginning, the 100% MR plateau is not affected until A/2 < M (see the curve with 
iV = 30, 50, and 100 in figure [6(a)). Then with further increasing N, the 100% MR 
plateau becomes narrow (see the curve with iV = 150 in figure E(a)), and it disappears 
at the width iV with its A = M. At last, when the width of ribbon is very wide with its 
A M, the MR is very small regardless of other parameters. This result is the same 
with recent work studying on the infinite wide graphene spin valve device [20]. So it is 
favorable to select the narrower nanoribbon to fabricate the spin valve. On the other 
hand, the length change of the central region does not affect the 100% MR plateau at 
all (see figure E^b)). Only the shapes of some subplateaus are slightly changed. 

In the above discussion, we assumed that the Dirac-point energy ec of the central 
region is equal to e a of the left and right leads. If ec departs from e Q , the conductance 
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Figure 6. (Color online) (a) The MR vs. Ep for different width N with the length 
L = 10 and (b) MR vs. E F for different length L with the width N = 50. The 
magnetization M = 0.01. 
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Figure 7. (Color online) The MR vs. Ep for different disorder strength W. The 
parameters are M = 0.02, N = 50 and L = 10. 



G in the antiparallel configuration still is zero when Ep G (-M, M) and < M < A/2, 
because the parity of transverse wave functions in the left and right FM leads can not 
match with each other. However, the conductance G for the parallel configuration is 
very large (> e 2 /h). So the 100% MR plateau is not affected by ec slightly departing 
from e a . 

How is the 100% MR plateau affected by the disorder? Here, we consider the 
Anderson disorder which exists only in the central graphene region. The on-site energy 
ec in the central region Hamiltonian becomes ec + Wi, where Wi is uniformly distributed 
in the range of [—W/2,W/2\. Figure [7] shows the MR vs. Ep at different disorder 
strength W, in which every MR curve for W ^ is averaged over up to 1000 random 
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configurations. The MR is reduced by the disorder, but the plateau shape is still hold 
and the value is about 40%-50% even in quite strong disorder strength. 




Figure 8. (Color online) The conductance G vs. the angle 9 for different 
magnetization M and different Fermi energy Ep. The size of the central region is 
N = 50 and L = 1 0. 

3.4- The conductance for the arbitrary angle 9 of the left and right magnetizations 

Now we analyze 9 dependence of the conductance. In figure [SJ the conductance G as a 
function of the angle 9 are plotted with combination of M and Ep. Here the shape of 
the curves G-9 is approximatively a cosine function regardless of the parameters M and 
E F . At E F = (i.e. at the center of the 100% MR plateau), the change quantity of the 
conductance G for the angle 9 from to 7r is about 2e 2 /h, which is quite large. About at 
9 = 2, the derivative dG/d9 reaches the extreme value, but at the region of < 9 < 1, 
the derivative dG/d9 is quite small. On the other hand, when Ep deviates away from 
the 100% MR plateau, the variation of the conductance G is also about 2e 2 /h or even 
larger for some values of Ep (see the curves for Ep = 0.05 in Fig.8), although the MR is 
about 50% there. However, at the region of MR = 0, the conductance G is independent 
of the angle 9 (not shown here). 

4. The case of armchair edge 

In this section, we investigate the MR in the armchair edge FM/normal/FM junction. 
In the case of the armchair edge, there are two types of energy band structure [25]. If 
the number of transverse atomic layers N is equal to 3m or 3m + 1 (m is an integer), 
an energy gap appears. This gap can be quite large for a narrow graphene nanoribbon. 
When the Fermi energy Ep is in the gap, the conductance is always very small regardless 
of the parallel or antiparallel configurations. So it is not interesting even if its MR is 
very large. For N = 3m + 2, the graphene nano-ribbon is metallic and its conductance 
is large. 
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Figure 9. (Color online) Left panel: The MR as a function of Ep in the case of 
the armchair edge with the width N = 41 and the length L = 40 atomic layers in the 
central region. Right panel: The energy band structure of an ideal armchair edge 
graphene ribbon with the width iV = 41. 

Figure [9] shows the MR as a function of Fermi energy Ep for the metallic armchair 
edge FM/graphene/FM junction with A 0.13i for N = 41. The conductance 
approaches 2e 2 /h for both the parallel and antiparallel configurations at a large range 
of Ep around zero energy. As a result the MR is very small regardless of M > A and 
M < A. On the other hand, when Ep is far away from (e.g. Ep > A), the MR can be 
over 50% in some regions of Ep. Anyway, the property of MR of the zigzag edge ribbon 
is much better than that of the armchair edge ribbon. So it is more favorable to fabricate 
the spin valve devices by using the zigzag edge ribbon. In experiment, the graphene 
nanoribbon with the specific parity edge has already been successfully fabricated [21]. 

5. Conclusions 

In summary, we have studied the electronic transport and magnetoresistance (MR) 
in graphene-based ferromagnetic/normal/ferromagnetic junction where the finite width 
and the arbitrarily relative orientation between the lead magnetizations are taken into 
account. For the zigzag edge case, the conductance for the parallel configuration 
is always larger than e 2 /h under any parameter condition, but for the antiparallel 
configuration the conductance is exactly zero when Fermi energy is around the Dirac 
point due to the band-selective rule. This leads to a 100% MR plateau. The 100% 
MR plateau is almost not affected by the length of the central graphene region. With 
the increasing of the width of the ribbon, the 100% MR plateau is kept stable in the 
beginning, but is smeared at the wide ribbon limit. In the presence of the disorder, the 
MR is slightly suppressed, but can still keep the plateau shape with the value about 50% 
even in quite large disorder strengths. When the orientation between magnetizations 
of two ferromagnetic leads is arbitrary, the conductance versus the relative angle 6 is 
similar to a cosine formation. What's more, for the armchair edge case, the MR is 
relatively small. Therefore, it is more favorable to fabricate the graphene-based spin 
valve devices by using the zigzag edge graphene ribbons. 
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